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Entangled photons are a crucial resource for quantum communication and linear optical quantum 
computation. Unfortunately, the applicability of many photon-based schemes is limited due to the 
stochastic character of the photon sources. Therefore, a worldwide effort has focused in overcoming 
^^ the limitation of probabilistic emission by generating two-photon entangled states conditioned on 

^^ the detection of auxiliary photons. Here we present the first heralded generation of photon states 

that are maximally entangled in polarization with linear optics and standard photon detection from 
? spontaneous parametric down-conversion 1 . We utilize the down-conversion state corresponding to 

the generation of three photon pairs, where the coincident detection of four auxiliary photons unam- 
biguously heralds the successful preparation of the entangled state ^ . This controlled generation 
of entangled photon states is a significant step towards the applicability of a linear optics quantum 
network, in particular for entanglement swapping, quantum teleportation, quantum cryptography 
and scalable approaches towards photonics-based quantum computing [3]. 
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Photons are generally accepted as the best candidate for quantum communication due to their lack of decoherence 
and their possibility of being easily manipulated. However, it has also been discovered that a scalable quantum 
computer can in principle be realized by using only single-photon sources, linear optical elements and single-photon 
^O detectors [4 . Several proof-of-principle demonstrations for linear optical quantum computing have been given, includ- 
t^^ ing controlled-NOT gates [MS], Grover's search algorithm [9l[T0], Deutsch-Josza algorithm [11^, Shor's factorization 
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algorithm [121 IE] and the promising model of the one-way quantum computation [14] . 
^D A main issue on the path of photonic quantum information processing is that the best current source for photonic 



entanglement, spontaneous parametric down- conversion (SPDC), is a process where the photons are created at random 

K^ times. All photons involved in a protocol need to be measured including a detection of the desired output state. 

/\ This impedes the applicability of many of the beautiful proof-of-principle experiments, especially when dealing with 

Cu multiple photon pairs [3 and standard detectors without photon number resolution. Other leading technologies in 

this effort are based on other physical systems including single trapped atoms and atomic ensembles [15], quantum 

dots [16 , or nitrogen- vacancy centers in diamond [17 . Although these systems are very promising candidates, each of 

these quantum state emitters faces significant challenges for realizing heralded entangled states; typically due to low 

coupling efficiencies, the uncertainty in emission time or the distinguishability in frequency of the photons created. 

However, the probabilistic nature originating from SPDC can be overcome by several approaches conditioned on 
the detection of auxiliary photons [2j [181 [E] • It was shown that the production of one heralded polarization-entangled 
photon pair using only conventional down-conversion sources, linear optical elements, and projective measurements 
requires at least three entangled pairs [20j . Here we describe an experimental realization for producing heralded two- 



photon entanglement along theses lines, suggested by Sliwa and Banaszek that relies on triple-pair emission from a 
single down-conversion source [2 . This scheme shows significant advantages compared to other schemes where either 
several SPDC sources and two-qubit logic gates [18] or more ancilla photons [19] are required. 

Current down-conversion experiments allow for the simultaneous generation of three photon pairs [2TH24] with 
typical detection count rates, dependent on the experimental configuration, of about 10~^ to 10~^ Hz. We use a 
setup of this kind such that the coincident detection of four auxiliary photons is used to predict the presence of two 
polarization-entangled photons in the output modes. The auxiliary photons thus herald the presence of a Bell state 
and it is not necessary to perform a measurement on that state to confirm its presence. 
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FIG. 1: Setup for the heralded generation of entangled photon pairs. Six photons are created simultaneously by exploiting 
higher-order emissions in a spontaneous parametric down- conversion process. The photons are passing a narrowband filter and 
are coupled to single-mode fibers. They are brought to beam splitters and the reflected modes are analyzed in \H/V) basis 
and in |±) basis, respectively, using polarizing beam splitters (PBS) and a half- wave plate (HWP) orientated at 45°. State 
characterization of the heralded photon pair in the transmitted modes is performed via polarization analysis and the help of 
quarter- wave plates (QWPs), HWPs and PBSs. 



Figure nigives a schematic diagram of our setup to generate the heralded state |0+) = 4= {\H)^^ \-^)t2 ^ l^)ti I^)t2)' 
where H and V denote horizontal and vertical polarization, respectively, whereas ti and ^2 correspond to the trans- 
mitted modes after the beam splitters. For generating the heralded state, |0^), three photon pairs have to be emitted 
simultaneously into spatial modes ai and a2, resulting in: 

|*3> = 1/2 • {\HHH)^^ \VVV)^^ - \HHV)^^ \VVH)^^ + \VVH)^^ \HHV)^^ - \VVV)^^ \HHH)J 

These photons are guided to non-polarizing beam splitters (BSi and BS2) with various splitting ratios. Our 
scheme only succeeds when four photons are reflected and measured in a four-fold coincidence. The two reflected 
photons of BSi are projected onto \H/V) basis for mode ri, while the two reflected photons of BS2 are measured in 
|±) = A={\H) ± |V)) basis for mode r2. We are interested in the case where one photon is present in each of the 
modes riH,iv and r2+,2-- Considering only these terms, the output state results in 

1*3) = c{e„e,) ■ \H)^^^ \v)^^^ 1+),^^ |-),^_ . i= {\H),^ \H),^ + \v),^ |y),J (1) 

where C(^i, 62) is a constant depending on the transmission coefficients of the beam splitters. The coincident detection 
of one and only one photon in the modes vih-, ^ly, r2+ and r2- heralds the presence of an entangled photon pair in 
state |(/)+) in the output modes ^1,^2- In the present scheme such a case can only be achieved by three-pair emission 
from SPDC. The contribution from two-pair emission is suppressed by destructive quantum interference in the HWP 
rotation used for r2+,2-. This quantum interference together with the use of number-resolving detectors ensures 
that the remaining two photons are found in the transmitted modes. If a high transmission of the beam splitters is 
chosen, it still can be assumed with high probability that the two photons are transmitted even without the use of 
number-resolving detectors. 

In our case of using standard detectors (PerkinElmer photo- avalanche diodes) the transmission of the non-polarizing 
beam splitters should ideally be as high as possible such that a measured four-photon coincidence corresponds to 
precisely four photons and thus heralds our desired state. Therefore for demonstrating this dependency we choose 
beam splitters with different transmission rates of 17 %, 50 % and 70 %. Obviously the disadvantage of increasing the 
probability of heralding a |(/)+) state - which in principle can be approximately unity - is a reduction in the four- fold 
coincidence rate for triggering this state. Only the recent improvements of laser sources enable stable UV beams 
with sufficient power for keeping the measurement time reasonable; in our case the typical counting time for one 
measurement setting varies from 24h to 72h. The actual rate R of the four-fold coincidences is about Ru/ss = 83, 
^50/50 = 14 and Rjq/^q = 0.4 per minute. 

For characterizing our heralded state, all polarization-dependent measurement outcomes in the output modes ti 
and ^2 that are triggered by the four- fold coincidence in modes vih^ ^ly, ^2+, and r2- are analyzed as a function 
of the beam splitter transmission. The measured probabilities of finding the various photon numbers in the output 
modes are shown in Table [l] and allow for the reconstruction of the diagonal elements of the density matrix in the 
photon number basis |ni)^ 1^2)^ , where rii and 7x2 are the Fock or photon- number states per spatial mode (see 
Supplementary Information). The dependency of the photon- number statistic on the beam splitter ratio can be 
clearly seen as the vacuum contribution Pq^o decreases with higher transmission rates. The resulting probability of 
heralded entanglement generation is graphically shown in Figure [2] as a function of the beam splitter transmission. 
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TABLE I: List of the photon- number probabiUties Pni;n2 of having m and n2 photons in the output modes ti and t2. The 
events with the same sum of photon numbers in modes ti and t2 are compared to the probabihty of obtaining one photon in 
each of the output modes. The error for each probabihty is fohowing a Poissonian distribution of the measured counts. See 
Supplementary Information for a more detailed analysis including polarization modes. 
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beam splitter transmission 

FIG. 2: Probability of heralded entanglement generation shown for various beam splitter transmissions. The deviation from the 
expected quadratic behavior (line) originates from high-order emissions, which increase the probability of measuring photon pairs 
in the output modes for higher beam splitter transmissions. The diamond-shaped data point originates from the experiment 
with the reduced laser-power and the error bars follow a Poissonian statistic. 
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These probabilities, defined as P = C6/(C4 • r]'^) where Cq (C4) is the six (four) -fold coincidence rate and r] is the total 
photon detection efficiency per mode, are Pu/ss = (2.5±0.2) %, P50/50 = (29.4±1.0) % and P70/30 = (77.2±6.6) % for 
the different transmission rates. Remarkably, these probabilities are achieved due to the high visibility of (86.2 ±0.7) % 
for the destructive four-photo quantum interference. 

The detection of more than one photon per spatial mode results from eight- or more-photon emissions due to the 
technical limitation of working in the high laser-power regime for optimizing count rates. Obviously, the detected 
six- fold coincidences for obtaining the Pi^i contribution also capture the higher-order emission, which cover about 
10 % of the coincidences for the high laser-power case. Naturally, the averaging of correlated and anti- or non- 
correlated measurement results decreases the quantum correlations of the desired output state. For fair and sensitive 
representation of this effect an additional set of polarization-correlation measurements is performed. Although these 
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beam splitter transmission 

FIG. 3: Experimentally obtained fidelities for the two-qubit polarization state with respect to the ideal state 0^) for various 
beam splitter transmissions (blue). The effect of higher-order emission is demonstrated by an additional experimental run 
with a reduced laser-power (green) . For this experiment a beam splitter transmission of 30 % is chosen for optimizing the 
required measurement time. As expected the probability of obtaining the heralded state 0^) increases with the transmittance 
of the beam splitters, whereas the polarization- state fidelities are not affected. The error bars are derived from Monte Carlo 
simulations based on a Poissonian distribution of the measured counts. 



polarization measurements are triggered by a four- fold coincidence, the requirement of obtaining an additional two- 
fold coincidence intrinsically leads to a post-selection of the two-photon polarization density matrix (see Figure |4|l and 
Methods). The corresponding fidelity, F^^^'^ of this measured photon pair with the corresponding entangled quantum 
state |^+) is Ff^^/g3 = (63.7 ± 4.9) %, F^^J^^ = (57.5 ± 3.4) %, and F^^J^^ = (61.9 ± 7.7) % for the different beam 
splitter ratios via local unitary transformations. 

For demonstrating the state fidelities' dependency on the laser-power (Figure [3|, an additional experimental run 
with a reduced laser-power of 620mW and beam splitter transmissions of 30 % was performed. The post-selected 
density matrix of this state in Figure ^p clearly shows an improvement of the polarization correlations, which is 
quantified by a fidelity of F^q%o ~ (^4.2 ± 8.5)%. From this data, we extract the tangle r [28 , a measure of 
entanglement that ranges from for separable states to 1 for maximally entangled states, as rso/70 = 0-55 ± 0.19. 
This density matrix, as commonly written in the coincidence basis, would allow a violation of local realistic theories 
by almost 2 standard deviations as it implies a maximum Clauser-Horne-Shimony-Holt [29l 130] Bell parameter of 
5* = 2.36 ± 0.22. This laser-power dependent noise is therefore not intrinsic in the setup and is only due to technical 
limitations, which can be overcome in future experiments by using photon-number discriminating detectors or with 
high-efficient down-conversion sources. 

These two-photon density matrices together with the measured photon number probabilities allow to calculate the 
state fidelity F'^^^^ = ((/)+| p |^^) of the output state including vacuum and higher-order terms. This measured total 
state fidelity can be extracted as F[^Jg^ = (0.0164±0.0010) %, F^J^^ = (0.0517±0.0029) % , F^^^^ = (0.176±0.013) % 
and F^J^Q = (0.497 ±0.041) %, which have significantly improved to standard down- conversion sources. These results 
suggest that with future gradual increases to the coincidence rate and the fidelity of entangled photons, the utilization 





FIG. 4: The effect of higher-order emission for the polarization correlations. The trade-off for the increased coincidence 
rates is manifested in the contribution of higher-order emission, a, Experimentally obtained polarization density matrix 
with a laser-power of 1.2 Watt and a beam splitter transmission of 50%. The captured eight-photon contribution leads to 
a background of a w~) state, b, The reduction of the background is demonstrated when reducing the laser-power. The 
experimentally reconstructed two-qubit polarization density matrix is measured with a laser-power of 0.62 Watt and a beam 
splitter transmission of 30%. The imaginary part of the density matrices is below 0.09 for all elements and hence not shown. 



of this scheme for quantum information processing tasks may not be far out of reach. 

This experiment presents the first feasible scheme for the generation of heralded entangled photon pairs with SPDC 
and single-photon detectors. This conditional method achieves a high preparation efficiency of up to 77% with 
measured fidelities of up to 84 % for the post-selected two-photon state. These results successfully underline its 
potential applicability for entanglement-based technologies. In conclusion, we highlight a multi-photon experiment 
that generates heralded entangled states as required for long-distance quantum communication and scalable quantum 
computing. We note that during the course of the work presented here we learned of a parallel experiment by 
Wagenknecht et al.[3T]. 
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and G. Mondl for assistance with the electronics. This work was supported by the Austrian Science Fund (FWF), 
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Commission under the Integrated Project Qubit Applications (QAP) and Quantum Interfaces, Sensors, and 
Communication based on Entanglement (Q-ESSENCE) and the 1ST directorate, the ERG Senior Grant (QIT4QAD) 
and the Marie-Curie research training network EMALI. 



METHODS 

We simultaneously produce six photons in the | ^3) -state by using higher-order emissions of a non-collinear type-II 
SPDC process. A mode- locked Mira HP Ti:Sa oscillator is pumped by a Coherent Inc. Verdi V-18 laser to reach output 
powers high enough to be able to exploit third-order SPDC emissions. The pulsed-laser output (r = 200 fs, A= 808 nm, 
76 MHz) is frequency-doubled using a 2mm-thick Lithium triborate (LBO) crystal, resulting in UV pulses of 1.2 W 
cw average. We achieve a stable source of UV-pulses by translating the LBO with a stepper motor to avoid optical 
damage to the anti-reflection coating of the crystal (count rate fluctuations less than 3% over 24 h). Afterwards, 
dichroic mirrors are used to separate the up-converted light from the infrared laser light. The UV beam is focused 
on a 2mm-thick /3-barium borate (BBO) crystal cut for non-collinear type-H parametric down- conversion. HWPs 
and additional BBO crystals compensate walk-off effects and allow the production of any Bell state. Narrowband 
interference filters (AA = 3 nm) are used to spatially and spectrally select the down-converted photons which are 
then coupled into single-mode fibers that guide them to the analyzer setup. There, the photon pairs are directed to 
non-polarizing beam splitters whose splitting ratios are chosen dependent on the experiment. The reflected modes 
then are analyzed in \H/V) basis and in |±) basis. At this specific angle, where the HWP rotates the polarization by 
45°, any possible four-photon state, emitted into the four modes, vih.iv and r2+,2-, will result only in a three-fold 
coincidence because of |H — )^^ = {\HH)^^ — \VV)^^)/^/2. Thus, these two photons will never be split up at the PBS 
and therefore never contribute to a fourfold coincidence detection. The typical photon coupling rates and detector 
efficiencies for each spatial mode are about 23% and 42 %. 

Our density matrix is reconstructed by a tomographic set of measurements, where combinations of the single photon 
projections \H/V)^ |±), and \R/L) = A={\H) ±i\V)), on each of the two photons in modes ti and ^2 are used. The 
most likely physical density matrix for our 2-qubit state is extracted using a maximum-likelihood reconstruction [25]- 
[27] . Uncertainties in quantities extracted from these density matrices are calculated using a Monte Carlo routine and 
assumed Poissonian errors. 
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SUPPLEMENTARY INFORMATION 

The analysis of the polarization-dependent photon- number distribution using the Fock basis is shown in Table |n| 
This data represents the diagonal elements of the density matrix, i.e. the probabilities of having n photons in each 
mode, starting from the vacuum contribution up to the case with one photon in each mode. 
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TABLE IL The Table lists the measured probabilities Pn^n ^n^v^'^2H ^n2v ^^ finding uih^ tiiv, ri2H and n2v photons numbers 
in the output modes tin Mv Mh and t2v with orthogonal polarizations H and V for the different beam splitter transmissions 
using our standard photo- avalanche diodes. The error for each probability is following a Poissonian distribution of the measured 
counts. 



For each beam splitter transmission ratio we used a (overcomplete) tomographic set of measurements for each 
basis (Jx^ cTy^ (Jz to reconstruct the polarization correlations of the (high-power) output state. In Table 
measured detection events that were triggered by a four- fold coincidence of the ancilla photons. 
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|0,1 


1,0) 


9 


13 


4 


15 


11 


7 


5 


11 


9 


|0,1 


0,1) 


5 


2 


6 


6 


2 


9 


5 


4 


10 


|0,0 


1,1) 


1 





1 


1 








2 


1 


1 


ll,l 


1,0) 





























ll,l 


0,1) 





























|1,0 


1,1) 





























|0,1 


1,1) 





























ll,l 


1,1) 





























30/70 


<Ti^V<^> 


CFy CFy 


ai^V?) 


ai^V?) 


ai^V^^) 


ai^)4^) 


ai^Vi^) 


4^Vi^) 


4^Vi^) 


|0,0 


0,0) 


53770 


59081 


61607 


54745 


55060 


57606 


26923 


57227 


59855 


|1,0 


0,0) 


568 


683 


677 


579 


584 


644 


257 


708 


752 


|0,1 


0,0) 


610 


571 


674 


605 


600 


587 


327 


587 


605 


|0,0 


1,0) 


506 


491 


551 


486 


450 


457 


220 


516 


570 


|0,0 


0,1) 


374 


463 


399 


380 


499 


442 


222 


408 


449 


ll,l 


0,0) 








1 


3 


2 


1 


3 


2 





|1,0 


1,0) 


19 


4 


17 


11 


11 


7 


5 


17 


13 


|1,0 


0,1) 


6 


10 


3 


5 


6 


14 


6 


10 


8 


|0,1 


1,0) 


10 


21 


4 


13 


6 


9 


3 


8 


5 


|0,1 


0,1) 


14 


3 


8 


2 


4 


8 


5 


7 


8 


|0,0 


1,1) 








1 


2 











1 


2 


ll,l 


1,0) 








1 




















ll,l 


0,1) 














2 














|1,0 


1,1) 





























|0,1 


1,1) 





























ll,l 


1,1) 






























11 



50/50 


ai^Vi^) 


ay ay 


ai'^ai'^ 


ai^V?) 




ai'^ai'^ 


ai^'ai'^ 


<yy CFx 


4^Vi^) 


|0,0 


0,0) 


32429 


15769 


16779 


14303 


14002 


16309 


31535 


17756 


15859 


|1,0 


0,0) 


956 


443 


390 


377 


387 


376 


838 


446 


405 


|0, 1 


0,0) 


749 


394 


482 


304 


292 


443 


884 


398 


374 


|0,0 


1,0) 


576 


277 


208 


204 


253 


314 


623 


338 


224 


|0,0 


0,1) 


627 


319 


433 


302 


292 


308 


720 


312 


381 


ll,l 


0,0) 


6 


1 


9 


5 


4 


5 


15 


1 





|1,0 


1,0) 


47 


3 


15 


13 


13 


16 


29 


20 


14 


|1,0 


0,1) 


25 


38 


9 


5 


14 


15 


22 


12 


8 


|0, 1 


1,0) 


19 


18 


8 


14 


4 


8 


21 


14 


7 


|0, 1 


0,1) 


33 


2 


11 


17 


9 


11 


48 


12 


11 


|0,0 


1,1) 


4 


1 


4 


3 








10 


2 


1 


ll,l 


1,0) 


2 








1 








2 








ll.l 


0,1) 





























|1,0 


1,1) 


1 

















1 








|0,1 


1,1) 





























ll,l 


1,1) 





























70/30 


ai^Vi^) 


ay ay 


ai^V?) 


ai^V?) 


(Jx CFy 


ai'^ai^^ 


ai'^ai'^ 


4^Vi^) 


4^Vi^) 


|0,0 


0,0) 


833 


1371 


1480 


1478 


1402 


3360 


3630 


1692 


972 


|1,0 


0,0) 


23 


55 


36 


56 


52 


94 


102 


64 


22 


|0, 1 


0,0) 


27 


43 


79 


55 


38 


136 


140 


57 


33 


|0,0 


1,0) 


26 


42 


48 


31 


45 


110 


121 


61 


25 


|0,0 


0,1) 


34 


47 


66 


52 


44 


147 


170 


76 


41 


ll.l 


0,0) 





1 


1 


1 


2 


2 











|1,0 


1,0) 


4 


1 


3 


3 


4 


6 


6 


4 


4 


|1,0 


0,1) 


2 


5 


4 


8 


4 


10 


3 


3 





|0, 1 


1,0) 


2 


4 


2 


1 


2 


11 


7 


3 


4 


|0,1 


0,1) 


4 





7 


1 


4 


12 


6 


4 


2 


|0,0 


1,1) 


1 








2 


2 








2 





ll,l 


1,0) 





























ll,l 


0,1) 














1 





1 








|1,0 


1,1) 





























|0, 1 


1,1) 





























ll.l 


1,1) 






























TABLE III: We show the measured photon-numbers for each polarization and spatial mode, Inii/, ni\^, n2i^, n2v), where ni/f, 
riiv denote the photon number for the orthogonal polarization states in output mode 1 and n^H^ n^v denote the orthogonal 



.(1)^(2) 
3 



by phase retarders in front of the polarizing beam splitters. 



